#---------------------------------------------------------------------
# To be copy-pasted in the initialization part of the RNPL source code
#---------------------------------------------------------------------

initialize sqrtming { [1:Nx][1:Ny][1:Nz] := (8/pi^3)* ( cos(pi*x/2) )^2 * ( cos(pi*y/2) )^2 * ( cos(pi*z/2) )^2 }

initialize gtt { [1:Nx][1:Ny][1:Nz] :=-1 }
initialize gtx { [1:Nx][1:Ny][1:Nz] := 0 }
initialize gty { [1:Nx][1:Ny][1:Nz] := 0 }
initialize gtz { [1:Nx][1:Ny][1:Nz] := 0 }

initialize gxx { [1:Nx][1:Ny][1:Nz] := (4/pi^2) * ( cos(pi*x/2) )^4 }
initialize gxy { [1:Nx][1:Ny][1:Nz] := 0 }
initialize gxz { [1:Nx][1:Ny][1:Nz] := 0 }

initialize gyy { [1:Nx][1:Ny][1:Nz] := (4/pi^2) * ( cos(pi*y/2) )^4 }
initialize gyz { [1:Nx][1:Ny][1:Nz] := 0 }

initialize gzz { [1:Nx][1:Ny][1:Nz] := (4/pi^2) * ( cos(pi*z/2) )^4 }

initialize FF  { [1:Nx][1:Ny][1:Nz] := 0 }
